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Abstract. We discuss the determination of fundamental parameters of 
'normal' hot, massive OB-type stars, namely temperatures, luminosi- 
ties, masses, gravities and surface abundances. We also present meth- 
ods used to derive properties of stellar winds - mass-loss rates and wind 
velocities from early-type stars. 
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1 Introduction 



IT) 

■ The vast majority of stars in the Milky Way and external galaxies have masses close 
\ to that of the Sun, and live long (>Gyr), relatively peaceful lives. In contrast, stars 

with initial masses in excess of «10M Q , with 0-type and early B-type spectral 
'"q! types on the main-sequence live for a relatively short life (^5-20 Myr), though 

■ have an astrophysical importance which belies their rarity. Massive stars and their 
O \ evolved descendants possess powerful stellar winds, such that they are important 

■ contributors to the chemical and dynamical evolution of galaxies, due to their 
\ influence on the composition and energetics of the interstellar medium (ISM). 

^ ■ Indeed, the primary origin of oxygen and other cv-elements in the universe is from 
massive stars. Ultimately, they end their life as a core-collapse Supernovae, causing 
further 'feedback' to the ISM. 



X 



Massive stars spend the majority of their lives in the blue part of the Hertzsprung- 
- - ■ Russell (HR) diagram, and so emit the majority of their energy in the far- and 
extreme- ultraviolet. Consequently, in normal galaxies they are the dominant 
source of Lyman continuum photons. A typical 07 dwarf emits 10 49 Lyman con- 
tinuum photons every second for its main sequence life, and so OB stars ionize 
material in their vicinity, producing ionized nebulae, or Hnregions. Within the 
Milky Way, the Orion Nebula is the closest birth place of massive stars. Further 
afield, 30 Doradus, the Tarantula Nebula represents the greatest single concentra- 
tion of massive stars within the Local Group of galaxies, containing many hundreds 
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of OB stars within its Giant H nrcgion. Further afield, starburst galaxies contain 
many so-called 'super-star clusters'. Massive stars are the only stellar objects to 
be observed in high redshift galaxies (e.g. MS 1512-cB58: Pettini et al. 2000) 

It is clear that massive stars play a major role in astrophysics. In this article I 
shall explain how one can determine the fundamental properties of main sequence 
massive stars. In the following article, I shall discuss post-main sequence massive 
stars and their end states. 

2 Basic Concepts 

The principal measured parameter of an individual star is its brightness in the 
night sky, and its colour, namely the difference between magnitudes at two or more 
wavelengths. Historically, photography was the method of choice for astronomers. 
Since photography is most sensitive to visible radiation, spectral classification 
traditionally concentrated on these wavelengths. Only recently have electronic 
detectors, such as CCDs supplanted photographic plates. 

Brightness, colour, and spectral types can lead to the luminosity and tempera- 
ture of the star, the fundamental parameters to place it on the Hcrtzsprung- Russell 
diagram. The apparent magnitude of a star is typically measured through a fil- 
ter, commonly the V defined in the Johnson UBV system. The absolute visual 
magnitude, My is related to the apparent magnitude, my via the relationship, 

My = m v - 5 + 5 log d - Ay 

where d is the distance (in parsec) and Ay is the U-band extinction, defined by 
Ay = R x E(B - V), where R ~3.1. E(B - V) is merely the difference between 
the observed colour excess mg — my = B — V in the Johnson system and the 
intrinsic (B — V)q colour. According to Wien's law, most of the stellar radiation for 
luminous hot stars with their ss 30, 000 K minimum T c g in emitted in the far UV 
regions. At optical wavelengths one is observing way out on the long wavelength 
tail of the spectral energy distribution (SED) for hot stars. In the UBV system, 
intrinsic colors of OB stars are thus nearly degenerate (e.g., differing only by a few 
hundredths of a magnitude), typically (B — V)o ~ —0.3 mag. 

The highly wavelength dependent extinction of starlight by interstellar dust 
dominates the colors of luminous stars in the Milky Way. The observed optical 
colors of hot stars are used primarily to determine properties of the intervening 
interstellar medium and are not very useful to evaluate stellar parameters. Us- 
ing the standard Galactic extinction curve of Seaton (1979), for every magnitude 
of extinction at the U-band, the far-UV extinction at 1200 A suffers four magni- 
tudes, whilst only 0.4 magnitudes is suffered at 1/nm and 0.1 magnitudes at 10/im. 
Consequently, regions, such as the Galactic Centre which are invisible at UV and 
optical wavelengths are reasonably transparent longward of the if-band. Inter- 
nal extinction in other galaxies differs somewhat from the Milky Way (Howarth 
1983; Bouchet et al. 1985). In addition to continuous dust absorption, there 
are also discrete interstellar absorption features, principally the atomic Lyman HI 
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and molecular H2 bands at UV wavelengths, whose strengths broadly scale with 
E(B - V) (see e.g. Pellerin et al. 2002) 

The m v measures the brightness of a star through a filter, but a luminosity 
is the emergent radiation over all wavelength bands. The Sun has an absolute 
visual magnitude of My=4.82 mag, whilst typical luminous hot stars have My — 
—5 mag, so they are intrinsically of order 10,000 times brighter in the visual. 
Considering their output over all wavelengths, they exceed this by a further factor 
of 10 to 100. In magnitudes, one can write, Mboi = My + B.C, where Mboi is the 
bolometric magnitude and B.C. is the bolometric correction. Thus latter quantity 
is strongly dependent on the T e g, having an appreciable value for hot, blue stars. 
One can relate the luminosity of a star in solar luminosities to Mboi via 

M bo] = -2.5 Iog(L/L ) + 4.74 

where L Q = 3.845 x 10 33 ergs -1 . 

In principal, one obtains the luminosity of a star by measuring its my and 
(B — V), obtaining a distance, adopting a B.C. and utilizing the above equations 
together with a calibration relating the spectral type of the star to its intrinsic 
B — V)o colour. A T c g will first need to be obtained so that the B.C. can be 
determined. Various B.C.(T e g) exist in the literature (e.g. Balona 1994). For O 
and early B stars, Vacca et al. (1996) derived 

B.C. = 27.66 - 6. 84 log T cff 

The dependence is clearly very steep for O stars, so small errors in T e g may lead 
to large errors in luminosity. For example, the stellar luminosity of an 03 dwarf 
with absolute magnitude My = —5 mag is 500,000L Q if T o ff=50,000K, but only 
370,000L Q , if T off =45,000K. 

It is, of course, no trivial matter to determine the distance of an individual 
hot star, and this represents the greatest difficulty with deriving luminosities for 
Galactic OB stars. Due to their large distances, direct determination via stellar 
parallax methods are only possible for a dozen cases with Hipparcos (Lamers et al. 
1997). Amongst the closest examples Hipparcos measured were ( Oph (09.5 Vn) 
and t Sco (BO. 2 IV) at distances of ~135 pc, a hundred times more distant from 
the Sun than the closest star, a Centauri. 

Within a decade, GAIA should enable reliable parallaxes to all optically visible 
Galactic OB stars. Until then, one has to resort to obtaining spectral type-My 
calibrations using OB stars in associations or clusters. Humphreys (1978) carried 
out a very extensive study for Galactic OB supergiants. This was extended to 
Magellanic Cloud associations by P. Massey and co-workers. These calibrations 
may then be applied to individual stars in the field. The advantage of Magel- 
lanic Cloud stars is that their distances are known to a precision better than 10%. 
Consequently, despite their relative proximity, luminosities of Galactic OB stars 
are more imprecise than those of stars in the Magellanic Clouds or beyond. The 
limitation of extragalactic stars is principally that of spatial resolution. Fortu- 
nately, one finds reasonably good consistency between calibrations from Galactic 
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and extra-galactic hot stars. Nevertheless, OB stars within a given spectral type 
show a substantial spread, which is thought to be typically ±0.5 mag. 

One may also use radio observations to study the Galactic distribution of mas- 
sive stars, since these are located in Giant Hnregions that are strong radio emitters 
(e.g. Smith ct al. 1978). OB stars are confined to the thin disk of the Milky Way, 
and (crudely) map out the spiral arm structure (Georgelin & Georgelin 1976). 
Knowledge of the ionizing fluxes of individual OB stars permits one to obtain to- 
tal numbers of stars in H urcgions. The ionizing fluxes of hot stars is a very strong 
function of temperature and to a lesser degree surface gravity. Vacca et al. (1996) 
provided a Teff calibration of Lyman (Qo) and Hei(Qi) ionizing fluxes for OB 
stars. The widely used T e g-spectral type relation for OB stars adopted therein 
(see below) has been substantially revised recently, such that Crowther & Dessart 
(1998) provide a calibration of Qo(2cff , log L), whilst Hubeny & Lanz (2003) have 
recently produced a comprehensive tabulation of ionizing fluxes for dwarf O stars 
at a range of mctallicities and gravities. 

3 Spectral Types 

The emergent continuum radiation from stars arises from the stellar photosphere 
while the absorption line spectrum comes from the atmosphere. In normal stars 
these regions overlap, can be treated as plane parallel, and are closely coupled to 
the stellar parameters. 

Numerical subtypes distinguish differences within the spectral letters and go 
from number to 9 (but with a few 0.5 divisions). Line widths and other line 
ratios are used to determine the atmospheric pressures (related to gas densities, 
local gravity, hence luminosity). These labels utilize roman numerals I-V, with the 
latter corresponding to the faintest luminosity class. While five luminosity classes 
are easily distinguished among cool K and M stars with their wide dispersion in 
My, only I, III, and V are adopted for O type stars which have a more limited 
My range. 

The most common system of spectral classification builds upon that of Morgan 
et al. (1943-MK) in which moderate spectral resolution is used to define natural 
groups of stars with similar spectral characteristics. This two dimensional system 
uses selected letters and numbers for the spectral types and subtypes, and roman 
numerals for the luminosity classes. Standard stars are typically selected for each 
spectral type and luminosity class. Classification of other stars is then carried out 
by obtaining their spectra and by using the spectral criteria an eye comparison to 
the standards. The dividing line for luminous hot stars at T e g just below 30000 K 
corresponds roughly to spectral type Bl on the Main Sequence. Thus O-early 
B spectral types are going to be the stars whose evolution we will follow in this 
article. 

The MK system for these stars was modified and defined by Walborn (1971a). 
This was based upon spectra taken in the blue spectral region, on photographic 
plates, with a (reciprocal) dispersion of 63 A/mm. His classification system ranged 
from 04 through B2.5. Walborn (1971b) next identified a group of early O stars 
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near 77 Carina, which extended his spectral subtype system in a natural way to 03. 
The classification for 0-early B type stars utilizes Hei and Hen line ratios, along 
with Sill, Sim and Siiv features. The leading lines are as follows: Hei: A4471 
and 4387 (all units here are A); Hen: A4542 and 4686; Sin: A4128-30; Sim: 
A4552; Siiv A4089 and 4116. Basically, the appearance of the A4542 Hen line 
distinguishes the hotter from the cooler B spectral types. The O subtypes are 
defined by the ratio A4471 He I/A4542 He 11 which monotonically decreases towards 
earlier spectral type, such that the A4471 line is very weak, or totally absent at 
03. In luminosity V class the A4686 Hen line behaves similarly to A4542 although 
it is a bit stronger and it persists into the earliest B type stars. 

For B spectral types, A4542 He II is absent and the He 1 lines are strong. Si iv 
lines, which are seen throughout the O types, begin to fade at BO. Sim lines 
appear and the Sim/Si,iv line ratio is useful between BO and B1.5. At B2 and 
later subtypes, Sin appears as Sim begins to fade and the Sin/Sin ratio comes 
into use. In B and late-0 type stars, the strength of these silicon features is 
luminosity dependent, being stronger in brighter stars. The helium lines, however, 
are less effected so that various silicon/helium line ratios are used as luminosity 
criteria. Among late-O-early-B stars, luminosity classes of la, lb, III and V are 
recognized. 

In most O type stars, though, a vastly different phenomenon begins to become 
apparent in the optical: the appearance of emission lines. Plaskett & Pearce (1931) 
noticed that in some O stars, the A4634-41 N III triplet came into emission along 
with A4686 He 11. These stars they labelled Of type. Walborn (1971a) pointed 
out that most O stars he observed showed the N III lines in emission even though 
A4686 remained in absorption or was not present. He labeled these as 0((f)) 
and 0(f) subtypes , respectively. We know now from detailed line modeling in 
O star atmospheres (e.g., Mihalas et al 1972) that Nm emission is a result of a 
fortuitous dielectronic recombination and implies nothing about the luminosity. 
On the other hand, the appearance of the A4686 He 11 line, whether in absorption, 
missing, or emission does say something about the stellar wind structure (Auer 
& Mihalas 1972) and the luminosity. The status of A4686 He 11 is useful as a 
luminosity criterion. Following Walborn (1971a) the luminosity class of all but 
the latest type O stars is as follows: if the line is in absorption, the type is O V; if 
missing, OIII; and if in emission, Of. The ((f)) and (f) nomenclature is (almost) 
redundant as nearly every O star, aside from the latest types, has A4640 N in in 
emission. 

Recently, the 02 spectral type has been added by Walborn ct al. (2002a), 
such that 02-4 stars are defined by the ratio of emission lines at A4058 N iv and 
A4634-41 Nm, in contrast with the usual He absorption line criteria for OB stars. 
The principal motivation for this was that all O stars for which 4471 He I was very 
weak or absent were previously grouped together as 03, yet exhibited a range of 
spectral morphologies, and stellar temperatures from quantitative analysis. The 
most massive, hottest stars known, with N iv 3> N ill emission, are re-defined as 
02 stars, and include HD93129A (02 If*) in the Carina Nebula. Until recently, 
this was thought to represent one of the most massive stars in the Milky Way 
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(Taresch et al. 1997) but it has recently been found to be a binary (Walborn 
2003). 

The MK and Walborn classification systems depend on the spectra of stars in 
the vicinity of the sun. These all have a composition similar to each other and 
to the sun itself. How would one classify stars that had a substantially different 
abundance? For O stars this is not a significant problem since the spectral criteria 
involve hydrogen and helium lines for the most part. Silicon ions do play a role 
in classification in the late O-early B stars and caution must be exercised in the 
cases of low metal abundance. For the LMC, with metal deficiencies of a factor 2 
or 3, the effects on the classification are barely noticeable. For the SMC, with an 
abundance down by a factor between 5 and 10, the silicon deficiency is obvious, 
but the silicon ion line ratios can be used with some confidence. Fortunately, the 
luminosity of these stars is known from their membership in the Magellanic Clouds 
and does not need to be established from their spectra. Lennon (1997) provides a 
revision to the spectral classification of B supergiants in the SMC. 

Walborn (1976) reviewed and clarified the observations that certain O-early 
B stars had anomalously strong nitrogen or carbon lines, but otherwise normal 
features. These he had labelled as OBN or OBC stars. The OBC stars were also 
notable for having weak nitrogen lines, but the OBN stars have more or less normal 
appearing carbon features. OBN stars are found among both main sequence stars 
and supergiants, while the OBC stars are exclusive to giants and supergiants. 
Other symbol or small letter suffixes have been introduced for anomalous OB 
spectra. These will be discussed in detail in the following article. A review of 
optical classification criteria for O-type stars is given by van der Hucht (1996). 

In addition to optical spectral classifications, it is also possible to assign, some- 
what cruder, ultraviolet or near-IR classifications. Detailed atlases of 1150-1800A 
ultraviolet OB star spectra taken with the IUE satellite have been compiled by N. 
Walborn and collaborators. These have recently been supplemented with far-UV 
datasets, covering the 900-1200A region with the Copernicus and FUSE satellites 
(Walborn et al. 2002b). Most recently of all, it has become possible to obtain clas- 
sification quality spectroscopy of OB stars in the near-IR, specifically the K-h&nA 
at 2^,m, which was carried out in a pioneering study by Hanson et al. (1996). This 
has proved to be an extremely valuable resource by which we can study individual 
hot stars in regions that are inaccessible to optical or UV observations. 

4 Stellar temperatures 

Stellar temperatures of early-type stars, essential for subsequent determinations 
of radii and luminosities, are derived from a comparison between observed pho- 
tometry or spectroscopy and models. 

4.1 Continuum techniques using LTE models 

The simplest case to consider is that for mid- to late- B dwarfs, with intermediate 
masses, which can be analysed using continuum energy distributions, via low dis- 
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persion spectroscopy or photometry. Stromgren uvbyfi photometry, coupled with 
Kurucz LTE line blanketed model atmospheres, provides a powerful technique for 
the determination of T e g in B-type and cooler stars via the Balmer discontinuity. 

Local Thermodynamic Equilibrium (LTE) Kurucz model atmospheres account 
very thoroughly for metal line blanketing and are widely employed for both early- 
and late- type stars. Although a detailed discussion of stellar atmospheres is 
beyond the scope of this book, line blanketing is the influence of thousands to 
millions of bound-bound spectral lines on the atmospheric structure. LTE means 
that the ionization state of the gas and the populations of the atomic levels can 
be obtained from the local T e and n e via the Saha-Boltzmann distribution, such 
that collisional processes occur faster than radiative processes. 

However, in the case of early- B dwarfs, B supergiants and O stars, several 
complications occur. Fatally for an approach based on the Balmer jump, this fea- 
ture disappears for stars with temperatures in excess of approximately 30,000K. 
Although the intrinsic UV-optical-IR spectral energy distribution docs differ be- 
tween BO and 02-3 stars, the effect is very subtle, and masked by uncertainties in 
interstellar reddening. In the case of B supergiants, the role of the Balmer jump 
is weakened, given that assumed solar abundances may no longer be appropriate. 
Consequently, the temperatures of such stars are derived via analysis of their line 
spectra, generally optical photospheric lines of helium (O-type) or silicon (B stars). 

4.2 Line techniques 

Generally, radiative processes generally dominate over collisional processes in hot 
star winds, so it is necessary to solve the equations of statistical equilibrium ev- 
erywhere, i.e. non-LTE. The problem with this is that a determination of popu- 
lations uses rates which are functions of the radiation field, itself is a function of 
the populations. Consequently, it is necessary to solve for the radiation field and 
populations simultaneously, which is computationally demanding, and requires 
an numerical iterative scheme to obtain consistency. Unfortunately, the problem 
is too complex for analytical solutions. Considerable effort has gone into devel- 
oping realistic non-LTE model atmospheres for early-type stars in recent years 
by a number of independent groups. Unfortunately, consistently treating metal 
'line blanketing' in extended non-LTE atmospheres is computationally demanding. 
Consequently, the earliest attempts by D. Mihalas in the late 1960's suffered from 
simplifications that only became surmountable two decades later, at which time 
large amounts of atomic data became available from the Opacity Project which 
was led by M. Seaton. 

From such techniques, temperatures of hot stars are obtained using ionization 
equilibrium techniques, i.e. from suitable diagnostic lines of adjacent ionization 
stages. Naturally, determinations of T e g from model atmospheres are critically 
dependent on assumptions going into those models and the accuracy with which 
surface gravities and elemental abundance ratios are determined. The majority of 
detailed analyses of OB stars have been subject to the assumption of hydrostatic 
equilibrium in a so-called 'plane-parallel' approach (e.g. Herrero et al. 1992; Puis 
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Fig. 1. Comparison between optical line profiles of HDE 269698 (04Iaf+) and synthetic 
spectra (dotted) calculated using an unblanketed, plane-parallel version of tlusty with 
T cff =46.5kK, logL/L =6.25, log 5 =3.7 (Crowther et al. 2002) 
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Fig. 2. Comparison between optical line profiles of HDE 269698 (04Iaf+) and syn- 
thetic spectra (dotted) calculated using the line blanketed, spherical code CMFGEN 
with T cff =40kK, logL/L Q =5.98, logg=3.6, M=8.5x 1O- 6 M yr" 1 , (3 = 1, Uoo=1750 
km s _1 (Crowther et al. 2002) 
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et al. 1996), where the spherical extension of the star is neglected. Optical Hei-n 
lines generally provide temperature diagnostics for O stars (e.g. Herrero et al. 
1992) while Sin-iv lines are used for B stars (e.g. McErlean et al., 1999). The 
typical range of temperatures spanned by O stars is T c g ~50,000K at the earliest 
subtypes to T off ~30,000K at BO V, and T cS ~10,000K at AO V. 

The temperature scale of OB dwarfs and supergiants has been compiled by 
various authors, most recently by Vacca et al. (1996). However, a word of cau- 
tion is necessary at this point. This, and all previous compilations are based on 
plane-parallel pure H-He models. The effect of line blanketing from millions of 
metallic lines on the emergent spectra and ionization balance was neglected in 
most literature papers. 

Within the past few years several non-LTE model atmosphere codes allowing 
for metal line blanketing have been developed, tlusty represents the current 
state-of-the-art for 'plane-parallel' studies of OB dwarfs (Hubeny & Lanz 2003), 
whilst spherically extended non-LTE line blanketed codes include CMFGEN (Hillier 
& Miller 1998), WM-BASIC (Pauldrach et al. 2001), fastwind (Santaloya-Rey 
et al. 1997) and the Potsdam code (Grafener et al 2002). Recent studies of OB 
dwarfs have been carried out using CMFGEN and TLUSTY, employing metal line 
blanketing, resulting in systematically lower temperatures (e.g. Martins et al. 
2002; Bouret et al. 2003). 

The case for O supergiants is further affected by the presence of strong stellar 
winds. These may contaminate ('fill-in') photospheric absorption lines, causing 
temperatures to be overestimated still further in the standard plane-parallel as- 
sumption (Schaerer & Schmutz 1994). Bohannan et al. (1990) demonstrated 
this for the O supergiant, £ Pup (04 If) in which wind blanketing reduced previ- 
ous determinations by >10%. Putting stellar winds and line blanketing together 
can cause downward revisions in T e g by up to 20% in extreme OB supergiants 
(Crowther et al. 2002; Herrero et al. 2002). This is illustrated in Figures 1 and 2, 
which show respectively, spectral fits to HDE 269698 (04Iaf+) in the LMC using 
a unblanketed plane-parallel code, revealing T e ff=46.5kK, in which the wind lines 
are neglected, and blanketed, spherically extended code, revealing T c s =40kK plus 
simultaneous fits to wind features in the optical, UV and far-UV (Crowther et 
al. 2002). Such major downward revisions have a strong influence on the stellar 
luminosity of individual stars, which of course scale as T^ s . 

5 Masses 

The mass of a star is perhaps its most fundamental property, yet it is generally 
very difficult to constrain for massive stars to better than a factor of two. This 
is because young massive stars of differing ages and masses cluster together very 
tightly in the upper left H-R diagram. There represent two, model dependent, 
methods that may be applied to derive a stellar mass, plus a third, more direct 
technique if the star happens to be a member of a double- lined, short period binary. 

In general, massive stars are not members of double-lined, short period binaries. 
The primary method of deriving a star's mass is via a determination of the surface 
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Fig. 3. Mass-luminosity relationship for young OB stars from binary orbits, taken from 
the Gies (2003) compilation, Hilditch et al. (1996) for DH Cep, plus Massey et al. (2002) 
for LMC stars 



gravity from fitting line profile wings of Balmcr hydrogen absorption lines. The 
spectroscopic mass is determined from this gravity determination, using 

GM , N 

In reality, an 'effective gravity' is determined for hot stars due to the complication 
of radiation pressure, as will be discussed later in this article, and is expressed as 
log<7 in cgs units. Typical values for OB dwarfs are \ogg =4-4.5, with somewhat 
lower gravities, implying more extended atmospheres, for supergiants, e.g. log g ^3 
in late O or early B supergiants. High dispersion optical spectroscopy is needed 
in order to model the Balmcr line wings accurately enough. 

Alternatively, the position of a star in the H-R diagram can be compared 
to predictions from evolutionary tracks to extract a stellar mass. This is most 
straightforward for main sequence stars, as it is often uncertain whether indi- 
vidual OB supergiants are evolving to cooler temperatures immediately after the 
main sequence, or blueward after an excursion to the red part of the H-R dia- 
gram. Surface abundances will generally allow one to discriminate between these 
alternatives, but these are difficult to reliably measure for hot stars. 

Herrero et al. (1992) first established a "mass discrepancy" between the sppec- 
troscopic mass determinations and those resulting from comparison with evolution- 
ary tracks. This remains to this day, with evolutionary masses up to a factor of 
two times higher than spectroscopic masses for OB stars. Potential problems lie in 
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both the stellar luminosities obtained from model analyses, see above, and in the 
evolutionary models. For the latter, there was previously a rather tight correlation 
between luminosity and mass for hot stars, but the recent inclusion of rotation into 
evolutionary calculations has revealed that alternative initial rotational velocities 
can lead to a considerable spread in (T e g, \ogL) for stars of identical initial mass. 

The most accurate method of determining masses for OB stars, and so estab- 
lishing which, if any, of these two approaches is correct, is if these are members 
of binary systems. Unfortunately, very few double-lined O+O systems have suf- 
ficiently short orbits to permit spectroscopic monitoring and determination of 
complete light curves practical. Gies (2003) presents a summary of recent stud- 
ies of detached (non-interacting) massive binaries. These permit a direct mass- 
luminosity relation for young OB stars, as shown in Fig. 3. Between 10 — 3OM , 
L oc M 2 - 5 , above which L oc M 1 ' 75 . 

Massey et al. (2002) have used very massive binaries in the LMC R136 clus- 
ter to confirm evolutionary mass determinations, including a mass of 57M Q for 
the 03 V primary of R136 #38 in close agreement with 53Mq from evolutionary 
models. Nevertheless, problems do remain in some cases. Consequently, the jury 
remains out for the present on this question. An added complication with the 
binary method is that stellar temperatures and bolometric corrections need to be 
derived or adopted from calibrations. Errors in such calibrations may exacerbate 
differences. Ultimately, the use of modern temperature calibrations in binary stud- 
ies and line blanketed stellar atmosphere models will likely resolve evolutionary 
and spectroscopic masses for early- type stars (Herrero et al. 2002). 

6 Stellar Wind Properties 

Winds are ubiquitous amongst massive stars, although the physical processes in- 
volved depend upon the location of the star within the H-R diagram. Mass-loss 
crucially affects the evolution and fate of a massive star, while the momentum and 
energy expelled contribute to the dynamics and energetics of the ISM. 

The existence of winds in O-type stars has been established since the 1960's, 
when the first rocket UV observations revealed the characteristic P Cygni signa- 
tures of mass-loss from Civ A1548-51, Siiv A1393-1402, and N v A1238-42 (Morton 
1967). A theoretical framework for mass- loss in hot stars was developed by Castor, 
Abbott & Klein (1975), known as CAK theory via line-driven radiation pressure. 
The fundamental characteristics of stellar winds are velocities and mass-loss rates. 
The former can be directly observed, whilst the latter relies of varying complexity 
of theoretical interpretation. 

6.1 Theory of radiatively driven winds 

The basic mechanism by which hot star winds are driven is the transfer of photo- 
spheric photon momentum to the stellar atmosphere through absorption by spec- 
tral lines. The velocity reached far from the star tends towards a terminal value, 
, due to the geometrical dilution of the photospheric radiation field, whilst the 
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mass- loss rate relates to the velocity field, v(r) and density, p(r) via the equation 
of continuity 

M = Airr 2 p(r)v(r) (6.1) 

for a spherical, stationary wind. 

The properties of stellar winds depend on both the number of metal lines 
available to absorb photon momentum, and on their ability to absorb, i.e. their 
optical thickness. Analytical solutions involve force multiplier parameters, a, S 
and k which characterize the outflow (Kudritzki et al. 1989). Of these, a controls 
the fraction of optically thick/thin lines, k the number of strong lines, and 5 the 
ionization balance. 



It can be shown that, 



1 — a 



(6.2) 



and 



M oc fcVCa-*) (6.3) 

From above, there is expected to be a close relationship between the terminal 
velocity and the escape velocity, 

v csc = {2g cff R)°- 5 (6.4) 

where the effective gravity 

GM . . . . 

5cff = -jp- (1 - T) (6.5) 

is the stellar gravity corrected for the reducing effect of radiation pressure (via 
Thompson scattering) on the gravitational potential via T which is related to the 
stellar luminosity and mass. It is clear that the determination of stellar escape 
velocities depends on knowledge of stellar masses and radii. 

Recent theoretical studies have indicated that CNO elements are principal line 
drivers for the outer, supersonic part of the wind, whilst iron group elements are 
responsible for the inner, subsonic part (Vink et al. 1999). The former determine 
the wind velocity, and the latter the mass-loss rate, at least for compositions close 
to Solar metallicity. 

Hot stellar winds are predicted to depend on the metal content, Z, as follows. 
Puis et al. (1996) showed that, k oc Z 1-< \ such that Ecm 16.31 implies 

M oc (6.6) 

Inserting typical values of a(=2/3) and suggests an exponent of ~0.6, 

close to 0.5 deduced by Kudritzki et al. (1987) and 0.7 derived by Vink et al. 
(2001). Since a is also predicted to be a function of metallicity, Leitherer et al. 
(1992) derived oc Z 013 for OB stars. 

Puis et al. (1996) and Kudritzki et al. (1999) obtained a theoretical rela- 
tionship between the bolometric luminosity of an early-type star and the wind 
momentum, Mv^ modified by the square-root of the stellar radius, i.e. the so- 
called wind-momentum relation, as follows: 
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M Voo {R/R Q ) 



0.5 



L l/(a-S) 



(6.7) 



For O stars, taking a = 0.6,6 — 0.05, implies a luminosity exponent of ~1.8, 
close to that observed for Galactic O supergiants (~1.5) and dwarfs (~1.6) accord- 
ing to Kudritzki & Puis (2000). The exponent for B supergiants is different, due 
to the change in ionization of the elements contributing to the radiative line accel- 
eration. Differences in coefficient follow from the recent revision in temperature 
calibration for O and early B stars (Herrero et al. 2002). This relationship pro- 
vides a potential independent method of determining distances to galaxies beyond 
the Local Group, although it is primarily of use for (less luminous) A supergiants 
which are visually the brightest hot stars in galaxies. 

6.2 Wind velocities 

Ultraviolet P Cygni profiles, ubiquitous in O-type stars provide a direct indication 
of stellar winds. Wind material approaching the observer within a column in front 
of the star is blucshiftcd by the Doppler effect. Scattering of radiation out of this 
direction causes a reduction in the observed flux on the blue side of the profile. 
From the regions on either side of this column wind material may scatter radiation 
towards the observer. This may occur from either the approaching (blue-shifted) or 
receding (red-shifted) hemisphere, the extra flux seen by the observer, originating 
from this scattered radiation, is a symmetric emission component on both sides 
of the profile. The overall effect is therefore, asymmetric with blue absorption 
and red emission. The wavelength of the blue edge of the absorption provides 
a measure of the asymptotic wind velocity. Accurate wind velocities of OB and 
W-R stars can be readily obtained in this way from HST or IUE observations of 
saturated Civ A1548-51, or Siiv A1393-1402 P Cygni profiles (Prinja et al. 1990; 
Lamersetal. 1995). The range in velocities encountered is considerable: Voo=3500 
kms -1 , corresponding to >1% of the speed of light, in the earliest O stars, to 100 
kins -1 in some AB supergiants. 

The number of OB stars observed in our Galaxy, principally with IUE, far 
exceeds that from external galaxies due to their relative UV brightness. Never- 
theless, sufficient extragalactic hot stars have now been observed, principally with 
HST, for us to be able to make useful comparisons. Wind velocities of LMC O 
stars differ little from Galactic counterparts (Garmany & Conti 1985), while a 
more prominent effect is observed in the SMC, particularly amongst early O stars 
(Walborn et al. 1995; Prinja & Crowther 1998). 

Theoretically, wind velocities of hot stars are expected to follow a so-called 
/3-law, 



where R is the stellar radius and (3 ~0.8-1.5. UV line profile modelling appears 
to confirm expectations. 




(6.8) 
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In spite of often uncertain stellar masses, Lamers et al. highlighted the so- 
called 'bi-stability' jump in wind properties around Bl (T c g ~21,000K), above 
which Woo — 2.65v esc , and « M ~ 1.4u osc . This is thought to result from the change 
in ionization of the dominant elements contributing to the line force (Vink et al. 
1999). 

6.3 Mass-loss rates 

Observationally, estimates of OB mass-loss rates are obtained from radio contin- 
uum fluxes, optical/IR line profiles, or UV wind lines. In extreme cases, M may 
exceed 10~ 5 Mq yr _1 in some early O supergiants, but values 10-100 times lower 
are more typical. 

6.3.1 IR radio continua 

Winds in hot stars can be readily observed at IR-mm-radio wavelengths via the 
free-free (Bremsstrahlung) 'thermal' excess caused by the stellar wind, under the 
assumption of homogeneity and spherical symmetry. Mass-loss rates can be de- 
termined via application of relatively simple analytical relations (e.g. Wright & 
Barlow 1975) which reveal that the continuum flux, S v oc v°- e , in a spherically sym- 
metric, isothermal envelope which is expanding at constant velocity. The emergent 
flux is set by the distance to the star d, mass-loss rate and terminal velocity as 
follows 



Other factors, including differences in composition, ionization balance, electron 
temperature generally play a minor role in normal OB stars. Barlow & Cohen 
(1977) used this approach together with IR excess fluxes to determine mass-loss 
rates for a large sample of Galactic OBA supergiants. Radio measurements of OB 
stars were carried out by Abbott and collaborators with the VLA, as summarised 
by Bieging et al. (1989), and by Leitherer and collaborators with the Australia 
Telescope (Leitherer et al. 1995). Unfortunately, OB stars with relatively weak 
winds do not show a strong IR excess or radio flux, so mass-loss results from this 
technique have solely been for nearby hot stars with dense winds. 

Collisions between stellar winds from stars in a binary system will cause non- 
thermal (synchrotron) radio emission, so care needs to be taken against overesti- 
mating mass-loss rates from (apparently) single stars in this way. Generally this 
is accomplished by carrying out observations at multiple radio frequencies. 

6.3.2 Optical and IR line profiles 

Ha has long been recognized as the prime source of mass-loss information in early- 
type stars. However, accurate determinations rely upon a complex treatment of 
the lines and continua, i.e. non-LTE, spherically extended models which treat the 




(6.9) 
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sub- and supersonic atmospheric structure. Such techniques have already been 
touched upon with regard to surface temperatures of hot stars. Problems with 
using fits to Ha to derive mass- loss rates may include blending with Hen A6560 
in O stars, together with nebular contamination from H n regions, and broadening 
of the central emission component by stellar rotation. 

Hot stars in the Magellanic Clouds provide us with the means to compare 
empirical mass-loss properties with predictions at lower Z . Although the number 
of stars analysed for the Magellanic Clouds remains small (Puis et al. 1996), the 
mass-loss rates of LMC and SMC O stars do appear to be lower than Galactic 
counterparts, in reasonable accord with theoretical predictions. 

Although moderate to high spectral resolution observations of OB stars are 
rare at near-IR wavelengths, there exist wind diagnostics analogous lines to optical 
lines (Bohannan & Crowther 1999). This provides the prospect of the determina- 
tion of mass-loss properties for hot stars obscured at optical wavelengths via IR 
diagnostics. 

6.3.3 UV P Cygni profiles 

UV P Cygni profiles from metal resonance lines provide direct evidence for mass- 
loss in hot stars. Unfortunately, mass-loss determinations require knowledge of 
elemental abundances, the degree of ionization of the ion producing the line, and 
knowledge about the form of the velocity law. 

Recent analyses of hot star winds via UV spectroscopy are generally carried 
out using the Sobolev with Exact Integration (SEI) method (Lamers et al. 1987; 
Haser et al. 1998), named after the pioneer of stellar wind studies, V. Sovolev from 
the late 1950's. In the SEI approach, UV P Cygni line profiles are fitted to reveal 
Mq, where q is the, a priori unknown, fractional population of a particular ion- 
ization stage within a particular clement. Direct ionization information is limited 
to observed ions from IUE or HST spectroscopy: Nv A1238-42, Siiv A1393-1402, 
and C iv A1548-51, so mass-loss rates may only be estimated via ionization balance 
information, generally from non-LTE atmospheric models. Far-UV FUSE obser- 
vations complement previous observations, providing information on many more 
ionization stages, including Nm A989-91, Cm A977, Siv A1062-72, Svi A933-44, 
Pv A1118-28. Massa et al. (2003) have used SEI modelling of line profiles mea- 
sured FUSE spectroscopy of LMC OB stars to derive empirical ionization balances, 
albeit relying on theoretical mass-loss rates from Vink et a. (2001). 

6.4 Structure 

There is extensive observational evidence that hot star winds are not smooth, 
steady outflows. Intensive monitoring of unsaturated UV P Cygni lines with IUE 
led to extensive variability in the absorption components. To produce such vari- 
ability, the wind structure must be on a large scale, covering a substantial fraction 
of the disk. In contrast, emission components of P Cygni profiles were found 
to be rather less variable, since it forms from a more global average of the wind. 
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Subsequent optical studies have also led to absorption and emission line variability. 

There is also indirect evidence that OB winds have a turbulent structure. Sat- 
urated P Cygni profiles have extended black troughs, which are thought to be a 
signature of a highly non-monotonic wind. Soft X-ray emission is also observed, 
which is thought to originate from embedded wind shocks (Chlebowski et al. 1989; 
Lucy & White 1980). Such shocks are also thought to be the means by which high 
ionization stages are observed, most notably Nv A 1238-42, in stars for which the 
T c ff is too low from normal means. This is via a process known as super- (or 
Auger) ionization. 

A likely explanation for this turbulent structure is theoretical evidence for a 
strong instability of line driving to small-scale velocity perturbations. There is a 
strong potential in line scattering to drive wind material with accelerations that 
greatly exceed the mean outward acceleration. Simulations demonstrate that this 
instability may lead naturally to a highly structured flow dominated by multiple 
shock compressions (e.g. Runacres & Owocki 2002). Such predictions are broadly 
consistent with recent Chandra X-ray line spectroscopy of £ Pup (Cassinelli et al. 
2001) and other early O stars. Interpretation of the X-ray spectroscopy of 8 1 C 
Ori indicates a magnetic origin. 

If hot star winds are structured, or clumped, it remains to be seen whether 
derived mass-loss rates of OB stars, as derived from optical or radio measurements, 
represent over-estimates of true values. Overall consistency between e.g. Ha and 
radio determinations (Lamers & Leitherer 1993) is rather good, suggesting that if 
winds are highly clumped, they are similar on scales of < 1.5i?* (Ha) and ^1000 
i?* (radio). Time-averaged line profiles of OB stars that have been intensively 
monitored do remain remarkably stable on scales of years. Nevertheless, recent 
FUSE spectroscopic observations of O supergiants in the Magellanic Clouds require 
that their winds are clumpy in order to simultaneously reproduce far-UV and 
optical wind lines, principally Pv A1118-28, unless phosphorus is undcrabundant 
relative to other elements in these galaxies (Crowther et al 2002). Evans et al. 
(2003) arrive at similar conclusions in late O and early B supergiants using Siv 
A1062-72. Although the ISM abundance of Phosphorus in the Magellanic Clouds is 
uncertain, Sulphur is well constrained. Consequently, the most credible possibility 
is that OB supergiants are clumped (see also Massa et al. 2003). 

7 Rotation 

During the formation process of all stars, including hot luminous OB stars, the 
conservation of angular momentum implies that they will commence their Main 
Sequence life rotating very fast, perhaps close to break-up velocity. Due to ro- 
tational broadening mechanisms, we are able to measure current OB rotational 
velocities from optical or UV photospheric line profiles. For OB stars, the princi- 
pal compilations are by Conti & Ebbets (1977) and Howarth et al. (1997). 

Due to projection effects, wsini is actually measured, where i is the inclina- 
tion of the system on the sky as viewed by the observer. In addition, measured 
projected rotation velocities underestimate the true projected equatorial velocity, 
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v e sini, due to limb darkening in rapid rotators (Collins & Truax 1995). From a 
sample of over 400 OB stars, Howarth et al. determined a median value of v sinz = 
90 kms -1 , which de- projects to a value 10% higher on average. The determination 
of rotational velocities is aided by some very slow rotators which serve as useful 
templates, r Sco (BO. 2 IV) is the most famous example, with vsini < 5km s _1 . 

One finds that the O Main Sequence distribution is much broader than that of 
supergiants, with only 3 supergiants amongst 33 stars with v r smi > 200 kms -1 . 
This can be understood in that as a star evolves from the ZAMS, its rotational 
velocity will decrease due to conservation of angular momentum, lowering the limit 
to the distribution. The various distributions are summarised by Howarth et al. 
(1997). 

Theoretically, Friend & Abbott (1986) first considered the effect of centrifugal 
acceleration on hot star winds, extended by Bjorkman & Cassinelli (1993) to allow 
for the azimuthal dependence, who found that M increases towards the equator, 
and Woo decreases, leading to an oblate 'wind-compressed disk'. This has been 
widely applied to Be stars, a subset of B stars exhibiting emission lines. Subse- 
quently, it was established by Owocki et al, (1996) that a quite different geometry 
may result, by additionally allowing for non-radial components of the line force, 
and the increased polar radiation flux via gravity darkening (von Zeipel 1924). 
Due to a reduced escape velocity, the equatorial velocity law is slower than at 
the pole, inhibiting the equatorial disk and the original scaling of mass-loss can 
be reversed, such that a prolate geometry results. By way of example, an early-B 
dwarf rotating at 85% of its critical, or break-up, velocity was found by Petrenz 
& Puis (2000) to have a strongly prolate structure (p po \ e / Pcq <15) with a polar 
(equatorial) terminal velocity of 1030 (730) kms" 1 . The global mass-loss rate 
was not found to deviate from its ID value by greater than 10-20%, except for 
supergiants close to the Eddington limit. 

8 Composition 

The French philosopher Auguste Comte famously, and inaccurately, predicted of 
stars in 1835 that we would never know how to study by any means their 
chemical composition .. ' Fortunately, theoretical knowledge of stellar atmospheres, 
atomic data, together with high-resolution spectrograph's now permits abundance 
determinations for most normal stars, including hot, luminous OB stars. 

This can only readily be carried out using optical photospheric lines of H, 
He and metals, based on the same non-LTE techniques described above that are 
used to derive T e ff, log g. Consequently, it is generally solely hydrogen- helium 
abundances which have been determined to date for a relatively large number of 
O stars (e.g. Herrero et al. 1992). Whilst most main sequence stars showed 
normal He abundances, He/H«0.1 by number, as expected, rapid rotators were 
surprisingly found to have enriched He at their surfaces, suggesting that rotational 
mixing has brought products of core nucleosynthesis to the surface while the star 
is still on the main sequence. Advances in model atmosphere modelling since 
then have generally supported these results. Recently, Bouret et al. (2003) have 
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analysed the CNO abundances of a sample of SMC O dwarfs using non-LTE line 
blanketed models. 

Studies of O and early-B supergiants reveal modest enrichment of He, whilst 
CNO elemental abundances have mostly been lacking until recently. Smith & 
Howarth (1994) studied late OC, O and ON supergiants, revealing increasing he- 
lium abundances for these, suggesting an evolutionary sequence in support of the 
morphological criteria. However, great care should be taken when analysing OB 
stars for abundance determinations since neglect of 'micro-turbulent' line broaden- 
ing may imply erroneous He-enrichments (Smith & Howarth 1998; McErlean et al. 
1998). There is clear evidence for enhanced nitrogen, accompanied with decreased 
carbon in OB supergiants from studies by Crowther et al. (2002), Hillier et al. 
(2003) and Evans et al. (2003). Notably, an ON9.7 Ia+ star in the LMC was 
found to be significantly more chemically processed than an otherwise identical 
09.7 Ia-I- star, whilst OC and BC stars were confirmed as CNO-'normaP. CNO 
abundances in OB supergiants in the SMC were found to exhibit similar degrees 
of chemical processing as A supergiants studied by Venn (1999). 

9 Further reading 

Kudritzki & Puis (2000) have provided an excellent Annual Review article on the 
properties of OB stars, whilst the book by Lamers & Cassinelli (1999) discusses 
Stellar Winds of both early- and late- type stars. Recent IAU Symposia involving 
massive stars include #212 (van der Hucht et al. 2003) and #215 (Maeder & 
Eenens 2003). 
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